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Abstract 

The Lovelock gravity consists of the dimensionally extended Euler densities. The 
geometry and horizon structure of black hole solutions could be quite complicated in this 
gravity, however, we find that some thermodynamic quantities of the black holes like the 
mass, Hawking temperature and entropy, have simple forms expressed in terms of horizon 
radius. The case with black hole horizon being a Ricci flat hypersurface is particularly 
simple. In that case the black holes are always thermodynamically stable with a positive 
heat capacity and their entropy still obeys the area formula, which is no longer valid for 
black holes with positive or negative constant curvature horizon hypersurface. In addition, 
for black holes in the gravity theory of Ricci scalar plus a 2n-dimensional Euler density 
with a positive coefficient, thermodynamically stable small black holes always exist in 
D = In + 1 dimensions, which are absent in the case without the Euler density term, 
while the thermodynamic properties of the black hole solutions with the Euler density 
term are qualitatively similar to those of black holes without the Euler density term as 
D > 2n + l. 
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Over the past years gravity with higher derivative curvature terms has received a lot 
of attention, in particular, in the brane world scenario and black hole thermodynamics. 
In the former, the motivation is to study the corrections of these higher derivative terms 
to the Newton law of gravity on the brane, or to avoid the singularity in bulk by using 
these higher derivative terms. In the field of black hole thermodynamics, it is expected 
to gain some insights into quantum gravity since the thermodynamic property of black 
hole is essentially a quantum feature of gravity. On the other hand, it is possible through 
the thermodynamics of black holes with AdS asymptotic to study the thermodynamic 
properties and phase structure of a ceratin field theory because due to the AdS/CFT 
correspondence, some higher derivative curvature terms can be regarded as the corrections 
of large N expansion of dual field theory. 

Among the higher derivative gravity theories, the so-called Lovelock gravity pQ is quite 
special, whose Lagrangian consists of the dimensionally extended Euler densities 

m 

£ = C ^n, (1) 

71=0 

where c n is an arbitrary constant and C n is the Euler density of a 2n-dimensional manifold: 
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Here the generalized delta function is totally antisymmetric in both sets of indices. Cq 
is set to one, the constant Cq is therefore proportional to the cosmological constant. Ci 
gives us the usual curvature scalar term. In order the Einstein's general relativity to be 
recovered in the low energy limit, the constant c\ must be positive. Here, for simplicity, we 
just set the constant c\ = \. The £2 term is the Gauss-Bonnet one, which often appears 
in the recent literature. Except for the advantage that the equations of motion of the 
Lovelock gravity, as the case of the Einstein's general relativity, do not contain terms with 
more than second derivatives of metric, the Lovelock gravity has been shown to be free 
of ghost when expanding on a flat space, evading any problems with unitarity [2]. Here 
it should be stressed that although the Lagrangian (0) consists of some higher derivative 
curvature terms, the Lovelock gravity is not essentially a higher derivative gravity since 
its equations of motion do not contain terms with more than second derivatives of metric. 
Just due to this, the Lovelock gravity is free of ghost [Sj. 

In the literature concerning on the Lovelock garvity, the extensively studied is the so- 
called Gauss-Bonnet gravity, whose Lagrangian is the sum of the curvature scalar term C\ 
and the Gauss-Bonnet term £ 2 > the Euler density of a 4-dimensional manifold. Sometimes, 
a cosmological constant is added to the Lagrangian. In this theory, the static, spherically 
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symmetric black hole solution was found in Refs. E|. The black hole solutions with 
nontrivial horizon topology were studied in Ref. Refs. [3 El discussed some aspects 
of holography of the Gauss-Bonnet gravity. In particular, it is worth mentioning here 
that with a positive Gauss-Bonnet coefficient C2, in spite of the asymptotic behavior 
(asymptotically (anti-)de Sitter 0IH] or flat P) of black hole solutions, it is found that a 
locally stable small black hole always appears when the spacetime dimension D = 5, which 
is absent in the case without the Gauss-Bonnet term, while D > 6, the thermodynamic 
behavior of the Gauss-Bonnet black hole is qualitatively similar to the case without the 
Gauss-Bonnet term (see also related discussions in [TTH ITT]). 

The Lagrangian looks complicated. It is therefore a little surprise to know that the 
static, spherically symmetric solution can be found in the sense that the metric function 
is determined by solving for a real root of a polynomial equation [3]. Since the gravity 

includes many arbitrary coefficients c n , it is not an easy matter to extract physical 
information from the solution. In Refs. [T2"l ITB"] by choosing a special set of coefficients, 
the metric function can be expressed in a simple form. These solutions could be explained 
as spherically symmetric black hole solutions. Black hole solutions with nontrivial horizon 
topology in this gravity with those special coefficients have also been studied in Refs. [TJJ 

EESj. 

For the general case with arbitrary coefficients c n , the static, spherically symmetric 
solution was found in Refs. jU E] 

ds 2 = -f(r)dt 2 + f(r)- l dr 2 + r 2 dVl 2 D , (3) 

where dQ 2 D denotes the line element of an (D — 2)-dimensional unit sphere and the metric 
function f(r) is given by 

f(r) = 1 -r 2 F(r). (4) 
F(r) is determined by solving for real roots of the following mth-order polynomial equation 

™ . . 16vrGM 

Here G is the Newton constant in D dimensions, Vlp = 2n^ D ^ 2 ^ 2 /T[(D — 2)/2] is the 
volume of an (D — 2)-dimensional unit sphere, M is an integration constant, and the 
coefficient c n is given by 

°° (D-l)(£>-2)' Cl ' 

c n = c n Tl 2 ™(D -i) for n > 1. (6) 
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The asymptotic behavior and causal structure of the solution have been analyzed in detail 
by Myers and Simon in Ref. [5]. It is found that even when the cosmological constant 
Co vanishes, the solution could be asymptotically de Sitter, flat, or anti-de Sitter, which 
depends on the coefficients c n ,(in other words, it can be seen from (J3J that the solution is 
asymptotically de Sitter, flat, or anti-de Sitter if F(r) approaches to a positive constant, 
zero or a negative constant as r — ► oo, respectively) and that (black hole/cosmo logical) 
horizon structure is quite rich. But we do not repeat them here. Further it is easy 
to conclude that the integration constant M is the ADM mass when the solution is 
asymptotically flat, while it corresponds to the AD mass for the asymptotically (anti-)de 
Sitter case [TT)j . 

Nowadays it is well-known that in asymptotically anti-de Sitter space, the black hole 
horizon could be topologically nontrivial: the horizon can be a closed hypersurface with 
positive, zero, or negative constant curvature [T2j. Such black holes are called topological 
black holes. Now we generalize the spherically symmetric solution (J3J) to more general 
case: 

ds 2 = -f(r)dt 2 + f(r)~ l dr 2 + r 2 dY? D , (7) 

where denotes a line element of an (D — 2)-dimensional hypersurface with constant 
scalar curvature (D — 2)(D — 3)k and volume Sx>. Here k is a constant. Without loss of 
generality, k can be set to ±1 or zero. In this case, the metric function /(r) becomes 

f(r) = k-r 2 F(r), (8) 

and F(r) still obeys the equation (jSJ) with Q D replaced by T, D . Following Myers and 
Simon [Oj, we can also make detailed analysis for the solution (JHJ) on its the horizon 
structure. Note that only when the solution is asymptotically anti-de Sitter, black hole 
horizon will appear for any fc; when the solution is asymptotical flat, it is possible to have 
black hole horizon only for the case k = 1; when the solution is asymptotically de Sitter, 
a cosmological horizon appears, of course, the black hole may also occur for k = 1 in this 
case. When the solution is asymptotically de Sitter, the solution (jHJ with any k is the 
generalization of the topological de Sitter spaces introduced in [T%] . 

For our purpose, without any detailed analysis, we just assume that a black hole 
horizon exists for the solution (J2J). Although the solution (j?j) looks involved, we will 
show that some thermodynamic quantities associated with the black hole can have simple 
expressions in terms of horizon radius. According to the metric (JZJ), the black hole horizon 
radius r + is determined via the equation /(r + ) = 0. Due to the equation (jHJ), one has 

r\ = k/F(r+). (9) 
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Note that when k — 0, it corresponds to r + = or F(r + ) = 0. For the former case, it 
indicates that the black hole horizon coincides with the singularity at r = 0. The Hawking 
temperature of the black hole can be obtained by using the periodicity of imaginary time 
in the metric. Continuing the black hole solution to its Euclidean section via r = it, 
the resulting manifold will have a conical singularity at the black hole horizon r + if the 
period (3 of the Euclidean time r is arbitrary. To remove the conical singularity, the period 
must be fixed to a special value. The periodicity of the Euclidean time appears in the 
quantum field's Euclidean propagator when one considers a certain quantum field in the 
black hole background. In quantum field theory at finite temperature, the period of the 
Euclidean time is explained as the inverse temperature, which is just the inverse Hawking 
temperature of black hole. For the black hole solution (J7J), the special value of the period 
of the Euclidean time is found to be 

P=l/T = 4n/f'(r)\ r=r+ , (10) 

where a prime denotes the derivative with respect to r. That is, the Hawking temperature 
is 

T = ^f(r)\ r=r+ = -^(2k/r + + r 2 + F' (r)\ r=r+ ) . (11) 



Here we have used the relation Q. To get a simplified expression, we have from the 
equation (JHJ) the mass of black hole in terms of the horizon radius 

lQnG ^ 

(12) 



When k = 0, it reduces to 



(D - 2)£ D r^ 1 , . 

To obtain F'(r)\ r=r+ , taking derivatives on both sides of equation (J5J) with respect to r 
and using (fT2"j) and we then get 



Substituting into (fTTj). we reach the expression of the Hawking temperature 



E n=o(D-^-^nk( kr- T »-■ 
4vrr+ En=inc n (kr- A 
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When k = 0, we find a very simple expression 

(D - l)c 

Tfc=0 = 4^ r+ ' (16) 

which is remarkable result: the Hawking temperature does not explicitly depend on other 
constants c n {n > 1). 

Another important thermodynamic quantity associated with black hole horizon is its 
entropy. Black hole behaves as a thermodynamic system, its thermodynamic quantities 
must obey the first law of thermodynamics, dM = TdS. Using this relation, in Ref. we 
have derived the black hole entropy in a higher derivative gravity theory, and in j^j we have 
obtained the same entropy of Gauss-Bonnet black holes as the Euclidean approach 
Here we use the first law to get the entropy of black hole (J7J). Integrating the first law, 
we have 

r r r + f)M 

s =J T ~ ldM =i r ~V r+ - (17) 

Here we have assumed that the entropy vanishes when the horizon radius shrinks to zero. 
Thus once the Hawking temperature and the mass are given in terms of the horizon radius, 
one can obtain the entropy of black hole using fTTj) . Substituting (JT2J) and (fTH|) into fTTj) . 
we arrive at 

Once again, the case with k = is very special. In that case, one can see from (fTHjl that 
there is only one term with n = 1 has the contribution to the entropy: 

S„ = (19) 

Note that C\ = 1 and S^r^ -2 is just the horizon area of black hole. We therefore conclude 
that in spite of the higher derivative terms, the entropy of black holes with k = always 
obeys the area formula of black hole entropy. For other cases with k — ±1, the area 
formula of black hole entropy does no longer hold obviously. 

Some remarks are in order here. First we notice that although the asymptotic behavior 
and horizon structure of the black hole solution ((7j) are complex, their thermodynamic 
quantities have simple expressions in terms of horizon radius. Their mass, Hawking 
temperature and entropy are given by (|T2|) . (fT5|) and (fT£|) . respectively. Second, when the 
horizon is a Ricci flat hypersurface, namely k = 0, the thermodynamic quantities of the 
black hole have quite simple forms given by (fTHj) . (fTBj) and (JTHJ), respectively. Further, 
from the relations (|T3|) and (jl6|) of mass and Hawking temperature to the cosmological 
constant Co, one has to have Co > 0, a negative cosmological constant, in order to make 
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these relations sense. In addition we see from (|T3|l and (JTHJ) that the black holes with 
k = are always thermodynamically stable with positive heat capacity. 

Third, as mentioned above, in spite of the asymptotic behavior of the Gauss-Bonnet 
black hole solution, the Gauss-Bonnet black holes with positive Gauss-Bonnet coefficient 
C2 always have a thermodynamically stable branch with small horizon radius in D = 5 
dimensions, while their thermodynamic properties are qualitatively similar to the case 
without the Gauss-Bonnet term if D > 6. Here we show that this feature persists for 
gravity with higher dimensional Euler density. For example, let us consider a gravity 
theory consisting of a cosmological constant term Co, a curvature scalar term R and a 
2n-dimensional Euler density C n . From (j!5|) . we have the Hawking temperature of the 
black hole 

_ (D - ljcgrg + (D- 3)r 2 + n - 2 + (D-2n- l)c n 

where we have already set k = 1. Suppose c n > 0, which makes the horizon radius r + have 
minimal value r + = (cf. jSl El El), we can easily see that the behavior of the Hawking 
temperature crucially depends on spacetime dimension. When D — 2n + 1, the Hawking 
temperature always increases monotonically from T = at r + = for small horizon 
radius, independent of the cosmological constant Cq. This is consistent with the fact that 
at much smaller scale than the cosmological radius 1 / yjfiol > (if h does not vanishes), 
the cosmological constant has a negligible effect on physics on that scale. Therefore in 
this case the small black hole is thermodynamically stable with positive heat capacity. 
Of course, for larger black holes, the behavior of Hawking temperature depends on the 
asymptotic behavior of the black hole solutions. From (|20|) it can be seen that the effect 
of the coefficient c n is small when r^ n ~ 2 > nc n . On the other hand, when D > In + 1, the 
Hawking temperature always decreases monotonically from infinity at r + = for small 
black holes, which implies that the heat capacity is negative, as the case without the Euler 
density term. For larger horizon radius, the effect of the Euler density term is once again 
small. Therefore the thermodynamic behavior of these black holes is qualitatively same 
as the case without the Euler density term as D > 2n + 1. 

Finally we see from the entropy (|T%|) that the first term is just quite familiar area term 
of black hole horizon, other terms comes from contributions of higher dimensional Euler 
densities. In this expression the cosmological constant term does not appear explicitly. 
This is an expected result since entropy of black hole is a function of horizon geometry j2Hj . 
Here we see that horizon topology also plays an important role for entropy of black holes 
in gravity with higher derivative terms. To see further the feature that black hole entropy 
is a character of horizon geometry and topology, let us add a Maxwell field to the Lovelock 
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gravity (|T|). We will see that entropy of the charged black hole in Lovelock gravity still 
have the expression ()18|). That is, the electric charge q does not appear explicitly in the 
black hole entropy expressed in terms of horizon radius. When a Maxwell field is present, 
we have a charged black hole solution with metric Here metric function f(r) is still 
given by ©, but F(r) has to satisfy PUT!] 

y a F n (r) - 16nGM q — (21) 

In this case, black hole horizon r + is still determined by the equation f(r + ) = 0. So the 
mass of black hole can be expressed in terms of horizon radius r + and charge q 

The Hawking temperature is found to be 

T Z™ =0 (D - 2n - ^kc^kr- 2 )- 1 - CD - 3)g 2 /rf - 6 
4vrr + E™=i^c n (A;r; 2 )«-i 

The variation of the mass (|22|) with respect to the horizon radius r + is 

(dM\ ( J D-2)S D r^- 4 ^ ^^-2xn-l , n o^/JIdV 



(23) 



9r i / IQnG 



J2(D-2n- l)kc n (kr- 2 ) n - 1 — (D — 3)g 2 /rf " 6 . (24) 

\n=0 ) 

Using (fTTj) once again, and keeping q as a constant in calculation, we get 

It has a same form as the entropy (fTS|) of a uncharged black hole. 

In summary we have first generalized the static, spherically symmetric black hole 
solution in Lovelock gravity to the case where black hole horizon can be a positive, zero or 
negative constant curvature hypersurface. Although the geometry and horizon structure 
of the black hole solution could be quite complicated, in terms of horizon radius, we have 
found that some thermodynamic quantities like the black hole mass, Hawking temperature 
and entropy, have simple expressions. In particular, the case with Ricci flat horizon is 
remarkably simple: these black holes are thermodynamically stable with a positive heat 
capacity and their entropy always obeys the horizon area formula. By explicit calculation, 
it has been shown that black hole entropy depends on not only the horizon geometry, but 
also the horizon topology structure. In addition, the feature has been found to be universal 
that for black hole solutions in gravity of Ricci scalar plus a 2n-dimensional dimensional 
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Euler density, when D = 2n + 1, the thermodynamically stable small black holes always 
appear with a positive heat capacity, which are absent in the case without the Euler 
density term. In D > 2n + 1, however, the thermodynamic properties of black holes 
become qualitatively similar to those of black holes without the Euler density term. This 
generalized the discussions of Gauss-Bonnet black holes to a more general case. 
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